y(4260) as a mixed charmonium-tetraquark state 
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Using the QCD sum rule approach we study the y(4260) state assuming that it can be described 
by a mixed charmonium-tetraquark current with J^^ = 1 quantum numbers. For the mixing 
angle around 6 ~ (53.0 ± 0.5)", we obtain a value for the mass which is in good agreement with the 
experimental mass of the y(4260). However, for the decay width we find the value Fy ~ (1.0 ± 0.2) 
MeV which is not compatible with the experimental value F ~ (88±23) MeV. Therefore, we conclude 
that, although we can explain the mass of the ^(4260), this state cannot be described as a mixed 
charmonium-tetraquark state since, with this assumption, we can not explain its decay width. 
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I. INTRODUCTION 

Many of the charmonium-like states recently observed 
in e+e" collisions by BaBar and Belle collaborations do 
not fit the quarkonia interpretation, and have stimulated 
an extensive discussion about exotic hadron configura- 
tions. The production mechanism, masses, decay widths, 
spin-parity assignments and decay modes of these states, 
called X,Y and Z states, have been discussed in some 
reviews [lH3)- Among these states, the F(4260) is par- 
ticularly interesting. It was first observed by BaBar col- 
laboration in the e~^e~ annihilation through initial state 
radiation and it was confirmed by CLEO and Belle 
collaborations 0. The r(4260) was also observed in 
the B- Y{4260)K- J/^Tr+n-R- decay Hi, and 
CLEO reported two additional decay channels: J/4'7r°7r° 
and J/^K+K- 0|. 

Since the mass of the y(4260) is higher than the 
Di*)D(*^ threshold, if it was a normal cc charmonium 
state, it should decay mainly to D^*^D^*\ However, the 
observed Y state do not match the peaks in e'^e~ — )• 
cross sections measured by Belle Q and 
BaBar [13, [HI- Besides, the *(35), ^{2D) and ^'(45') 
cc states have been assigned to the well established 
^'(4040), \E'(4160), and ^'(4415) mesons respectively, 
and the prediction from quark models for the 'if {3D) state 
is 4.52 GeV. Therefore, the mass of the ^(4260) is not 
consistent with any of the 1 cc states @, S[l2|- 

There are many theoretical interpretations for the 
y(4260): tetraquark state [l^, hadronic molecule of 

DoD* [3, XciUiM, XciP M, J/^foi^) [13. a 
hybrid charmonium |18| , a charm baryonium |l9| , a cusp 
[2(3l - |23 |. etc. Within the available experimental informa- 
tion, none of these suggestions can be completely ruled 
out. However, there are some calculations, within the 
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QCD sum rules (QCDSR) approach (sl. [23l425t . that can 
not explain the mass of the y(4260) supposing it to be a 
tetraquark state [2^, or a DiD, DqD* hadronic molecule 
[Hi, or a J/^/>/o(980) molecular state 

In this work we use again the QCDSR approach to the 
y(4260) state including a new possibility: the mixing 
between two and four-quark states. This will be imple- 
mented following the prescription suggested in [2^ for 
the light scalar mesons. The mixing is done at the level 
of the currents and was extended to the charm sector 
in Ref. [2^, in order to study the X(3872) as a mixed 
charmonium-molecular state. In particular, in Ref. (29| . 
the mass and the decay width of the X(3872), into 27r 
and Stt, were evaluated with good agreement with the 
experimental values. Agreement with the experimental 
results has been also obtained, applying this same ap- 



proach, in the study of the X(3872) radiative decay I30|, 
and also in the X{3872) production rate in B decay [31 1. 

In the next sections we consider a mixed charmonium- 
tetraquark current and use the QCDSR method to study 
both, mass and decay width, of the y(4260). 



II. CONSTRUCTING THE TWO-QUARK AND 
FOUR-QUARK OPERATOR 

In order to define a mixed charmonium-tetraquark cur- 
rent we have to define the currents associated with char- 
monium and four-quarks (tetraquark) states. For the 
charmonium part we use the conventional vector current: 

J'm^^ = Ca{xh^,Ca{x), (1) 

while the tetraquark part is interpolated by (26| 

{Qlix)Cj5Ct{x)){qd{x)jf,'y5CcJ{x)) -\- 



'(4) ^abc^de 



V2 



(2;) C75 7m Cfe (a;) ) (gd (a;)75 Cc J (x) ) 



(2) 



As in Refs. |28l. |29|. we define the normalized two-quark 
current as 

1 



(3) 
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and from these two currents we build the following mixed 
charmonium-tetraquark J^'-^ = 1 current for the 
y(4260) state: 

j,{x) = Me) ]'^\x) + com jj^H^), (4) 



III. THE TWO-POINT CORRELATION 
FUNCTION 



To obtain the mass of a hadronic state using the 
QCDSR approach, the starting point is the two-point 
correlation function 

U^Aq) = ^ I d^x e''--(0| r[j^ (x)jt (0)] |0) 



Iloiq'f-^, (5) 



where jfi{x) is the mixed charmonium-tetraquark inter- 
polating current defined in Eq. (U). The functions I{i{q^) 
and no((7^) are two independent invariant functions as- 
sociated with spin-1 and spin-0 mesons, respectively. 

According to the principle of duality, Eq. (O can be 
evaluated in two ways: in the OPE side, we calculate the 
correlation function in terms of quarks and gluon fields 
using the Wilson's operator product expansion. The phe- 
nonienological side is evaluated by inserting, in Eq. ([SJ, 
a complete set of intermediate states with 1 quantum 
numbers. In this side, we parametrize the coupling of 
the vector state Y with the current defined in Eq. Q 
through the coupling parameter Xy 



{Q\j,{x)\Y)=XYe^. 

where is the polarization vector. Using Eq. ([( 
can write the phenomenological side of Eq. ([S]) as 



id) 



M2 



(6) 



we 



(7) 



where my is the mass of the Y state and the dots, in the 
second term in the RHS of Eq. ([7]), denotes the higher 
resonance contributions which will be parametrized, as 
usual, through introduction of the continuum threshold 
parameter sq (s^ . 

The OPE side can be written in terms of a dispersion 
relation 



n 



OPE 



id' 



„OPE 



is) 



(8) 



where p'~'^^ (s) is given by the imaginary part of the cor- 
"~~(s) = Im[UOPE(s)]. In this 



relation function: 



np 



OPE 



side, we work at leading order in as in the operators and 
we consider the contributions from the condensates up to 
dimension eight. Although we will consider only a part of 



the of the dimension 8 condensates (related to the quark 
condensate times the mixed condensate), in Ref. |33| it 
was shown that this is the most important dimension 8 
condensate contribution. 

Considering the current in Eq. ^ , Eq. ^ in the OPE 
side can be written as 



^ cos' {6) Ulliq) + szn' (9) UtUq) 



M sin (9) cos {9) 
v2 



TilUq) + nfM 



(9) 



with 



^^J e''-(0|T[j;(x)j^T(o)]|o). (10) 

Clearly Il^j^{q) and Il'Jj^{q) are, respectively, the correla- 
tion functions of the J/tp and [eg] [eg] tetraquark state. 

After making a Borel transform in both sides, and 
transferring the continuum contributions to the OPE 
side, the sum rule in the structure for the vector 
meson can be written as 

A|,e-"-/*^« = M!cos2(0) nf (M|) + s^n2(0) nf (A/|) 



2 

M 

V2 



sin{9)cos{9) 



nf(Afi) + nr(Afi) 



(11) 



where 



nf(M^) = 



ds c-^/''^hllM 



n 
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(G2 



4mJ 



(12) 



nf(M|)= J dsc-^/''^-UUs)+pfl^^{s)+ 

Ami ^ 
r44 

1(8) 



(13) 



Til\Ml)= j dse-^/''--pll^{s)+nll^^^{M\ 

im'i 



(14) 



The expressions for the spectral density p{s) appearing 
in Eqs. (jl2p - (|T^ for the charmonium and tetraquark 
states, as well as the mixed terms are listed in Appendix 
A. 

By taking the derivative of Eq. pTjl with respect to 
l/Mg and dividing the result by Eq. (fTT|) . we obtain 



m 



dK{Mg 6) 
2 _ <i(l/M^) 



(15) 
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where 



K{Ml,0) = M!cos2(0)nf (M|) + sm2(0)nf (Ml) 



2 



^sin{e)cos{e) 



nf(M|) + nf(Af|) 



Eq. (|T5|) will be used to extract the mass of the 
charmonium-tetraquark state. 



A. Numerical Analysis 

In Table U we list the values of the quark masses and 
condensates that we have used in our numerical analy- 
sis. For a consistent comparison with results obtained 
for the others works using QCD sum rules, these param- 
eters values used here are the same values used in Refs. 



TABLE I. Quark masses and condensates values. 



Parameters 


Values 


mc{mc) 


(1.23 ± 0.05) GeV 




-(0.23 ±0.03)^ GeV^ 


{qga.Gq) 




ml 


(0.8 ±0.1) GeV^ 




(0.88 ± 0.25) GeV' 



The continuum threshold is a physical parameter that, 
in the QCDSR approach, should be related to the first 
excited state with the same quantum numbers. In some 
known cases, like the p and J/V', the first excited state 
has a mass approximately 0.5 GeV above the ground 
state mass. Since in our study we do not know the ex- 
perimental spectrum for the hadrons studied, we will fix 
the continuum threshold range starting with the smaller 
value which provides a valid Borel window, as explained 
below. Using this criterion, we obtain sq in the range 
4.6 < < 4.8 GeV. 

Reliable results can be extracted from the sum rule 
if is possible to determine a valid Borel Window. Such 
Borel window is obtained by imposing a good OPE con- 
vergence, the dominance of the pole contribution and a 
good Borel stability. To determine the minimum value of 
the Borel mass we adopt the criterion for which the con- 
tribution of the higher dimension condensate should be 
smaller than 15% of the total contribution. Thus, M|^j„ 
is such that 



OPE summed up to dim n-l(A/|,j 



total contribution(il/|„^-jJ 



= 0.85. 



(16) 



In Fig. [T] we plot the relative contributions of all the 
terms in the OPE side. We have used -Jsq = 4.70 GeV 
and 9 = 53°. For others 9 values outside the range 



52.5° < 9 < 53.5°, we do not have a good OPE con- 
vergence. From this figure we see that the contribution 
of the dimension-8 condensates is smaller than 15% of the 
total contribution for values of M| > 2.4 GeV^, indicat- 
ing a good OPE convergence. Therefore, we fix the lower 
value of A/| in the sum rule window as: M|^j„ = 2.4 
GeY\ 




2.5 3.0 3.5 

Ml (GeV2) 



FIG. 1. The OPE convergence in the region 2.0 < M| < 6.0 
GeV^ for = 4.70 GeV. We plot the relative contribu- 

tions start with perturbative contribution (line with circles), 
and each other lines represents the relative contribution after 
adding of one extra condensate in expansion: +{qq) (dot- 
dashed line), +{G'^) (long-dashed line), +{qga.Gq) (dotted 
line), +{qq)^ (dashed line) and {qq){qga.Gq) (solid line). 



To determine the maximum value of the Borel mass 
i^'^Bmax) must analyse the pole-continuum contribu- 
tion. Unlike the pole contribution, the continuum con- 
tribution increases with due to the dominance of 
the perturbative contribution. Therefore, the maximum 
value of the Borel mass is determined in the point that 
the pole contribution is equal to the continuum contri- 
bution. 

In Fig. [2] we see a comparison between the pole 
and continuum contributions. It is clear that the pole 
contribution is equal to the continuum contribution for 
Ml = 2.90 GeV^. Therefore, for ^ = 4.70 GeV^ and 
9 = 53° the Borel window is: 2.4 < Af| < 2.90 GeV^. 

After we have determined the Borel window, we can 
calculate the ground state mass, which is shown, as a 
function of M|, in the Fig. [3l From this figure we see 
that there is a very good stability in the ground state 
mass in the determined Borel Window, which are repre- 
sented, through the crosses in Fig. [21 

Varying the value of the continuum threshold in the 
range = 4.70 ± 0.10 GeV, the mixing angle in the 
range 9 = (53.0 ± 0.5)°, and the other parameters as 
indicated in Table I, we get: 



my 



(4.26 ±0.13) GeV, 



(17) 
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FIG. 2. The pole contribution (divided by the total, pole 
plus continuum, contribution) represented by solid line and 
the continuum contribution (dotted line) for the ^/sq = 4.70 
GeV. 
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FIG. 3. The mass as a function of the sum rule parameter M|, 
for ^ = 4.60 GeV (dotted line), ^ = 4.70 GeV (sohd line), 
y/so = 4.80 GeV (long-dashed line). The crosses indicate the 
valid Borel Window. 



IV. THE VERTEX FUNCTION AND THE 
DECAY WIDTH OF THE y(4260) 

The QCDSR technique can also be used to extract cou- 
pling constants and form factors. In particular, in Ref. 
|36| the authors determined the form factors and coupling 
constants in many hadronic vertices containing charmed 
mesons, in the framework of QCD sum rules. In this sec- 
tion, we will use the QCDSR approach to determine the 
coupling constant associated with the vertex YJ/il^a to 
estimate the decay width of the process Y — >■ J/ipTnr. We 
are assuming that the two pious in the final state come 
from the a meson. 

To determine the coupling constant associated with 
the vertex YJ/tpa^ we must evaluate the vertex function 
(three-point function) defined as 

n^,(p,p',g)= j d^xd^ye'P'--e"i-yU^,{x,y), (19) 

with p = p' + q and n^,y(a;, y) given by 

n,,ix,y) = {0\T{f;{x)r{y)j^H0)m. (20) 

The interpolating fields appearing in Eq. ((20)) are the 
currents for J/ip, a and y(4260), respectively. The cur- 
rents for J/tp and Y were defined by Eqs. ([1]) and 
For the meson a, we have 



r = -^{Ua{x)Ua{x) + da{x)da{x) 



(21) 



As in the case of two-point function studied in the pre- 
vious section, the three-point correlation function defined 
by Eq. can also be described in terms of hadronic 

degrees of freedom (Phenomenological side) or in terms 
of quarks and gluons fields (OPE side). In order to eval- 
uate the phenomenological side of the sum rule we insert, 
in Eg. p^ . intermediate states for Y, J/ip and a. Using 
the definitions: 

{0\jt\J/i'ip'))=m.^Ue^ip'), 



\,r\a{q))^A,, 



which is in a very good agreement with the experimental 
mass of the ^(4260). 

Once we have determined the mass, we can use this 
value in Eq. pT|) to estimate the meson-current coupling 
parameter, defined in Eq. We have used the same 

values of the sq, 9 and Borel Window used for the mass 
calculation. Thus, we get: 



Ay = (2.00 ± 0.23) X 10"^ GeV^ 



(18) 



The parameter Xy gives a measure of the strength of the 
coupling between the current and the state. The result in 
Eq. (fT8)) has the same order of magnitude as the coupling 
obtained for the X{3872) [U, for example. 



(y(p)ijrio) = Aye:(p), 

we obtain the following relation: 



n|fj'^"'(p,p',g) 



(p2 _ to|,)(p'^ - TO^)(g2 - ml) 



(22) 



(23) 



where the dots stand for the contribution of all possible 
excited states. The form factor, gYtpa(<f')^ is defined by 
the generalization of the on- mass-shell matrix element, 
(J/?/'cr|F), for an off-shell a meson: 

{J/il)a\Y) = gY^a{q^){p ■P^{p)-<^{p)-p' ■<^{jp)p-(*{jp')), 

(24) 
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which can be extracted from the effective Lagrangian that where Vaj3 = daYp — dpYa and A"^ = d^ip^ — d^ip", are 
describes the coupHng between two vector mesons and the tensor fields of the Y and iJj fields respectively, 
one scalar meson: 

C = tgri^aV^fsA"^ a (25) 
I 




FIG. 4. Diagrams which contribute to the OPE side of the sum rule for the structure p'^q^. 



In the OPE side, we work at leading order in as and we 
consider the condensates up to dimension five, as shown 
in Fig. m We have chosen to work in the p'^g^^ structure 
since it has more terms contributing for the OPE. Taking 
the limit = p''^ = —P^ and doing the Borel transform 
to — >■ M^, we get the following expression for the sum 
rule in the structure p'^g^: 



(26) 

where = —q^, and B(Q^) gives the contribution to the 
pole-continuum transitions (29ll37l-l39t . n^o^^' (M^, Q^) 
is given by 



n(o^^)(Af2,g2) = 



sin{0) 



mdqgcr.Gq) rl - 2a(I - a)i {glC^) 



a{\ — a) 



25 



(27) 



The sine present in Eq. (|27p indicates that only the 
tetraquark part of current in Eq. ^ contributes to the 
OPE side. In fact, the charmonium part of the current 
gives only disconnected diagrams that are not considered. 



In Eq. (|26ll m^, and f^, are the mass and decay con- 
stant of the J and m„ is the mass of the tr meson. 
Their values are: = 3.1 GeV, = 0.405 GeV [13, 
and rria = 0.478 GeV [4l|. The parameters Ay and 
represent, respectively, the coupling of the Y and a states 
with the currents defined in Eq. ^ and (|22p . The value 
of Ay is given in Eq. (fT5)) . while A^, was determined in 
Ref. HI and its value is A„ = 0.197 GeV^. 

Similarly to what was done to get my in Eq. (|15p , one 
can use Eq. ([^5]) and its derivative with respect to 
to eliminate B{Q^) from these equations and to isolate 
SY'ipa-iQ'^)- A good sum rule must be as much indepen- 
dent of the Borel mass as possible. Therefore, we have to 
determine a region in the Borel mass where the form fac- 
tor is independent of M^. In Fig. [5] we show gy^jaiQ^) as 
a function of both and . Notice that in the region 
7.0 < < 10.0 GeV^, the form factor is clearly stable, 
as a function of AP ^ for all values of . 

The squares in Fig. [6] show the dependence of 
gYi,a{Q ), obtained for AP = 8.0 GeV^. For other val- 
ues of the Borel mass, in the range 7.0 < AP < 10.0 
GeV^, the results are equivalent. Since we are interested 
in the coupling constant, which is defined as value of the 
form factor at the meson pole: = ^^a' '^^ need to 
extrapolate the form factor for a region of where the 
QCDSR is not valid. This extrapolation can be done by 
parametrizing the QCDSR results for gvipa-iQ'^) using a 
monopole form: 



6 



9i 

92 + Q'' 



(28) 



We do the fit for ^ = 4.74 GeV. We notice that 
the results do not depend much on this parameter. The 
results are: 

51 = (0.58 ± 0.04) GeV; g2 = (4.71 ± 0.06) GeV^. (29) 




FIG. 5. gYtpaiQ^) values obtained by varying both and 



The solid line in Fig. [S]sho"ws that the parametrization 
given by Eq. (pS)) reproduces very well the QCDSR re- 
sults for gviiaiQ^), in the interval 2.0 < < 4.0 GeV^ 
where the QCDSR is valid. 

The coupling constant, gYiixr is given by using = 
-ml in Eq. We get: 




FIG. 6. QCDSR results for gviiaiQ'^), as a function of Q^, 
for ^/so = 4.76 GeV (squares). The solid line gives the 
parametrization of the QCDSR results through Eq. 



The decay width for the process y(4260) — > J/i/'cr - 
J/tpTTTT in the narrow width approximation is given by 



^{Y J/i,nn) = 
as 87rm 



2ml- 



rc,(s)mo 



pis) 



TT {s - ml)^ + {mcyT^{s))^'' 
with p(s) given by 



(31) 



Pis) = ^ 5 —, (32) 

Zmy 

where A(a, &, c) = + 6^ + — 2ab — 2ac — 26c, and (s) 
is the s-dependent width of an off-shell a meson [4l| : 



5y0a = gy^ai-ml) = (0.13 ± 0.01) GeV^^ 



(30) 



The error in the coupling constant given above comes 
from variations in sq in the range 4.6 < sq < 4.8 GeV^, 
and in the mixing angle 52.5° < 9 < 53.5°. 

In Table II, we show the other values of the coupling 
constant corresponding to the values of y/so that we have 
considered in our calculations. 

Table II: Monopole parametrization of the QCDSR results 
for the chosen structure, for different values of ^/so 



(GeV) gy^^(Q') (GeV-^) gy^ajQ'' = -^ral) (GeV"^) 



4.6 
4.7 
4.8 



Q^+4.66 

0.57 
Q^+4.71 

0.53 
0^+4.77 



0.14 
0.13 
0.12 



T.is) = Tj ';^T-r\ ^^ (33) 
y A{mp,mj^,mj^j s 

where Fqct is the experimental value for the decay of the a 
meson into two pions. Its value is Tqc ~ (0.324 ± 0.042 ± 
0.021) GeV 

The invariant amplitude squared can be obtained from 
the matrix element in Eq. ([231) • We get: 



l-^P = 9Yi,ais)fimY,m^,s), 



(34) 



where gy^ais) is the form factor in the vertex YJ/'ipa, 
given in Eq. (P5|) using s = —Q^, and 

/(my, m^, ■5) = (^Yf^l + ^(m^ + ml - s)A . 
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Therefore, the decay width for the process y(4260) 
J l^-KTT is given by 



r = J, 



(35) 



where we have defined 

1=1 ds g^^„{s)T„{s){w?Y~m^^ + s) 

J(2m„)2 

p{s) 



X f{mY,m^,s)- 



(s-m2)2 + (m,r,(s))2 



(36) 



Hence, taking variations on sq and 9 in the same intervals 
given above, we obtain from Eqs. (j30p - ([55|) the following 
value for the decay width 



r(r ^ J/i/mTT) = (1.0 ± 0.2) MeV, 



(37) 



which is not compatible with the experimental decay 
width value expected for the F(4260) state which is 



around F. 



± 23) MeV 



Ilf^{Mg) state as well as the mixed terms Ilf^{M^) and 
n|2(Mg). The contributions for the last two are equal, 

that is, nf(M|) = nf (M|). 

For the charmonium contribution, the spectral densi- 
ties are written below 12411 



pert / \ 

P22 (s) 



237r2 



nr(Ml).-i^gM! 



(1 + 2m\l s)^J\^^^^^^^^ 



(Al) 



da.{ 2 



m2(l - 7q - 2q2) 4m^ 
a(\-afMl ^ M4(l-a)3/ 



e -B 



(A2) 



For the tetraquark we have 



SUMMARY AND CONCLUSIONS 



In summary, we have used the QCDSR approach to 
study the two-point and three-point functions of the 
K(4260) state, by considering a mixed charmonium- 
tetraquark current. In the determination of the mass, 
we work with the two-point function at leading order in 
as and we consider the contributions from the conden- 
sates up to dimension eight. A very good agreement with 
the experimental value of the mass of the F(4260) is ob- 
tained for the mixing angle around B w (53.0 ± 0.5)". 

To evaluate the width of the decay y(4260) — > J /il^nn, 
we work with the three-point function also at leading 
order in and we consider the contributions from the 
condensates up to dimension five. We assume that the 
two pions in the final state come from a cr meson. The 
obtained value for width is Fy « (1.0 ± 0.2) MeV, which 
is much smaller than the experimental value: T^xp ~ 
(88 ± 23) MeV. 

Therefore, we conclude that the y(4260) exotic state 
cannot be described as a mixed charmonium-tetraquark 
state. 
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Appendix A: The spectral densities for charmonium 
and tetraquark 

Next, we list all the spectral densities that appear 
in Eqs. ([r^ - fH]) for charmonium n22(Mg), tetraquark 



pert/ \ 



da 
a^ 



d/3 



F^{1 ~a 



/3)2 - 3F(1 + a + /3) j , (A3) 

(A4) 



P44 



(S) = - 



32 • 2" 



da 



d(3 



2m^a(l -a- Pf - 3mlF{l - a - (3) 
{2a^ + a(8 + 3/3) + ^(1 + 13) - 2) + 
-H6F2/3(1 -2a- 213) 



(A5) 



l-a 



-a(l + /^) - 2/3^ 



da 



16mi + 2H 



1 — a 



^(^ml{9-3a-5l3) + 7F 



(A6) 
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32 • 2%2 

(qq) {qGq) 
3 • 2V 

ml{qq){qGq) 
3 • 247r2 

q;2 — 2m2 
M|q;(1 - a) 



(1 - m-nll s)^l - Ami/ >iM) 



daa(5-6a) (A8) 



da X 



(A9) 



Finally, for the mixed term we have 



3 • 237r2 



(1 + 2m2/s)^l -4m2/s, 



(AlO) 



n<f'>(M^) = - 



^2) ^ ml{qq){qGq) f da 



3-237r2 _y 





(All) 



In all these expressions we have used the following def- 
initions: 



F ={a + I3)ml - aj3s, 



and the integration limits are: 





Ami/ s 


2 




Am^/s 


2 






am^ 



[sa — 77l2) 
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